Singular fate of the universe in modified theories of gravity 
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In this paper we study the final fate of the universe in modified theories of gravity. As compared 
with general relativistic formulations, in these scenarios the Friedmann equation has additional terms 
which are relevant for low density epochs. We analyze the sort of future singularities to be found 
under the usual assumption the expanding Universe is solely filled with a pressureless component. 
We report our results using two schemes: one concerned with the behavior of curvature scalars, and 
a more refined one linked to observers. Some examples with a very solid theoretical motivation and 
some others with a more phenomenological nature are used for illustration. 
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I. INTRODUCTION 

Refined astronomical observations of luminosity dis- 
tances derived from Type la supernovae provide reliable 
evidence of the current cosmic speed up of the Universe 
(see Q for the pioneering results and for the latest). 
In fact, such measurements are the only direct indication 
of that phenomenon (see for e.g. Q), but at the same 
time they are complementary with other key observations 
such as those of the CMB spectrum and the global matter 
distribution. Explaining this surprising behavior in the 
large-scale evolution of the Universe represents a major 
theoretical problem in cosmology, and several approaches 
have been coined to try and provide a compelling answer 
to this riddle. 

The main stream approach is to consider the Uni- 
verse is filled with an exotic fluid, known as dark energy 
S 0, but then one also has to demand the cosmic 
soup (made of dark energy and the rest of components) 
has some Goldilocks properties to comply with the obser- 
vations. Alternatively, the idea that cosmic acceleration 
might be due to modifications to general relativity has 
received considerable attention as well (see [13] for re- 
views and for specific modifications). In such frameworks 
models displaying cosmic acceleration could be devised 
with less fine-tuning and unnaturality as compared to 
general relativistic dark energy scenarios Specula- 
tions in the direction of modified gravity are, in princi- 
ple, legitimate as there are no cosmological tests probing 
scales as large as the Hubble radius. We only have reason- 
able evidence of the validity of the gravitational inverse 
square law up to 300 Mpc (through the ISW effect) p^ . 
However, the Hubble radius is two orders of magnitude 
larger, so our large-scale tests on general relativity are 
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not stringent enough. 

The additional degrees of freedom of these various set- 
tings, as compared to the standard picture of cosmology 
prior to the revolution ignited in 1998, have given rise to 
a collection of new cosmological evolutions with bold fea- 
tures, future singularites being the most perplexing ones. 
In this respect, attempts to classify somehow the sort of 
future singularities to be expected in new devised cosmo- 
logical evolutions are of interest. A popular classification 
route in the literature [l^ relies exclusively on proper- 
ties of the curvature tensor and scalar quantities derived 
from it. From that perspective, a number of new terms 
in cosmology, such as the celebrated "big rip" [l^j have 
been coined to designate extremality events associated 
with blow-ups of scalars constructed from the curvature 
tensor, along with less popular ones like "quiescent sin- 
gularities" "sudden singularities" "big brake" 
[19| or "big freeze" [l^l (the number of names is larger 
than the actual name of different extremality events). 

Now, even though treatments of singularities in the 
fashion of 1^ are of interest, there are subtle and most 
relevant properties inherent to cosmic evolution which 
can only be unveiled through the more sophisticated 
consideration of observers (see [13, [lH for a detailed 
account). Indeed, curvature is a static concept, as it 
is only provides information of what happens at each 
event. Conversely, information retrieved from tracking 
the observers along their trajectories is more dynamical 
in nature, and therefore more enlightening if carefully 
analyzed. Interestingly, this scheme allows discussing 
whether the singularities encountered are weak or strong. 
Thus, if one's ultimate goal is to draw rigorous con- 
clusions about the final fate in the Universe, both ap- 
proaches are, in our view, complementary. 

In this paper we address the problem of future sin- 
gularities in modified gravity cosmologies. We examine 
carefully the interrelation between the modifications and 
the singularities to be expected, and we try and give a 
unified vision by reporting our results using the scheme 
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concerned with the behaviour of curvature scalars [13| 
and the one grounded on observers [11] ■ 

Ideally, modifications of general relativity should be 
derivable of a parent theory allowing for a covariant for- 
mulation of full-fledged field equations, otherwise, neither 
density perturbations nor solar system predictions could 
be computed. This is, actually, an aspect of the problem 
which does not affect our discussion, as we only work at 
the level of the Friedmann and energy conservation equa- 
tions. Whenever the literature offers relevant examples 
for which the underlying theory is known, we will use 
them to illustrate our findings, but, occasionally, we will 
also resort to phenomenological examples. 

The plan of the paper is as follows. We propose a 
perturbative formulation of the Friedmann equations, for 
which two cases are distinguished depending on whether 
there is a critical energy density (which affects the form of 
the formulation). Then, we calculate the corresponding 
asymptotic expression of the scale factor, and bulding on 
earlier works we present our classification. We round up 
the dissertation with relevant examples and summarize 
in the last section. 



II. MODIFICATIONS OF FRIEDMANN 
EQUATION 

There have been many attempts to modify Einstein's 
theory of gravity from different points of view in order to 
cope with the observed acceleration of the expansion of 
the universe. One possibility arises from modifications to 
the Einstein-Hilbert action leading to the so called f{R) 
gravity theories (many aspect of this theoretical setup 
have been recently reviewed in plj |. The equations gov- 
erning the large-scale geometry of the Universe in such 
settings are of fourth order in the metric approach, and, 
on top of that, for f{R) gravity theories to evade compati- 
bility issues with observational tests complicated models 
are required (22| (see however [23| for a different per- 
spective). Mild applications of Ockham's razor princi- 
ple, combined in graceful cases with physical motivations, 
have lead to the consideration that contending modified 
gravity schemes could perhaps be more advisable. This 
is the case of the proposals originated by assuming the 
Universe is a 3-brane embedded in a higher dimensional 
bulk. 

Instead of grounding our discussion in specific theoret- 
ical frameworks, we propose a perturbative expression for 
the Friedmann equation of an expanding universe, which 
intends to comprise most of the models in the literature. 

With this aim in mind, we write a modified Friedmann 
equation in the form 




= H^ = hoip - p,)^" +hi{p- p,)«i + • • • . (1) 



Thus, we assume the squared Hubble factor can be ex- 
pressed as a power series in the density p of the matter 
content of the universe around a specific value p, , for 
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TABLE I: Singularities in FLRW cosmological models 



which a qualitative change of behavior is expected. The 
exponents are real and ordered, < S,i < ■ ■ ■ ■ The 
coefficient /ig is obviously positive. 

The equation system is closed by assuming, in addi- 
tion, the validity of the usual energy conservation equa- 
tion 

p + 3H{p + p)=0. (2) 

The perturbative formulation represented by Eq. ([T]) 
can accommodate the Friedmann equations of the ex- 
isting modified gravity proposals with a known parent 
covariant theory, as well as others with a phenomenolog- 
ical origin. Note, as well, that the A-cold dark matter 
(LCDM) or cosmic concordance scenario [l^l is trivially 
comprised within this framework: 

= ho + hi{p- p,), 

with ^0 = 0, ^1 = 1, hi = SttG/S, and /iq — ^iP* = A/3, 
so, actually, the parameter is not fixed. 

The main purpose of the modifications is to provide an 
accelerated evolution of the Universe without resorting to 
an exotic fluid, so it is usually assumed the Universe is 
simply filled with cold dark matter (p = 0) , and this will 
be our working hypothesis as well. In this case, the en- 
ergy conservation equation ^ can be straightforwardly 
integrated: 

^ = -3- ^ pa^ = K, (3) 
p a 

which gives a one-to-one map between the energy density 
and the scale factor through the integration constant K. 

If we perform a power expansion of the scale factor in 
time, 

a(t) = co|t-ior° +ci|<-iir'i + ••• , 

where the exponents rji are real and ordered, rjo < rji < 
• • • , following [l^ we shall be able to classify the sin- 
gularities encountered at a time to. It is expected that 
most models allow this sort of expansion. However, there 
are models arising in loop quantum cosmology [25| which 
show accelerated oscillations that fall out of this scheme, 
though most of our conclusions can be extended to them. 

The classification of singularities in weak and strong 
follows the ideas of Ellis and Schmidt : The curvature 
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may be finite or infinite at one event, but what is physi- 
cally relevant is whether free-falling (or even accelerated) 
observers meet the singularity in finite proper time [13] ■ 
It is clear that if they take infinite time in reaching the 
curvature singularity, this would be indetectable. 

Furthermore, if instead of ideal uncxtcndcd observers 
we consider finite objects, the key issue is whether tidal 
forces at the singularity are strong enough to destroy 
them or weak, so that there could be objects that would 
survive beyond the singularity. This would mean that the 
weak singularity is by no means the end of the universe. 

Following this ideas, Tipler [2^ modelled extended ob- 
jects by three perpendicular vorticity-free Jacobi fields 
travelling along a causal geodesic and forming and or- 
thonormal frame with the velocity u of this. If the 
geodesic hits a singularity in finite time and the volume 
spanned by a set of three three such vectors remains fi- 
nite, the singularity is considered weak, since an object 
is not crushed. Tipler argues that in this case the met- 
ric could be generically extended beyond the singularity. 
Otherwise, if the volume is not finite for every set of vec- 
tors, the singularity is considered strong. 

Thinking of cosmic censorship conjectures, Krolak [2^ 
suggested and alternative definition of strong curvature 
singularities that just required diminishing, instead of 
vanishing, volume of the finite object and is therefore 
easier to comply. 

Compliance with these definitions for FLRW models 
can be checked resorting to integrals of the Ricci tensor 
along causal geodesies with respect to proper time r [soj : 

If this integral diverges at a value tq 



[ dr' [ dT"Riju'u\ 
Jq Jq 



(4) 



the causal geodesic meets a strong curvature singularity 
according to Tipler's definition. 

And for Krolak's definition divergence of this other in- 
tegral at To 



Jq 



(5) 



means meeting a strong curvature singularity. 

The application of these results to FLRW models is 
summarized in Table I, which is a simplified version of 
the one in |,14i] . 

The last columns refers to the classification of future 
singularities in fTs^l): 

• Type I: "Big Rip" : divergent a, p, p. 

• Type II: "Sudden" : finite a, p, divergent p. 

• Type III: "Big Freeze" : finite a, divergent p, p. 

• Type IV: "Big Brake" : finite a, p, p, but divergent 
higher derivatives. 



Even though the modifications will only lead to ac- 
celeration for certain values of the parameter our 
forthcoming discussion on the asymptotic behavior of the 
scale factor is valid for any value of f^nd it only relies 
on the ordering of the exponents. For this reason, our 
scheme comprises as well modifications to gravity which 
are not able to explaining the current acceleration, such 
as, for instance, the non-self-accelerating branch of the 
DGP modification scenario. 

Inserting the modified Friedmann equation into the 
conservation equation one gets: 



P 
P 



3 hi 



2 x/hc 



(6) 



Our purpose it to integrate the latter by considering all 
the possibilities which arise from different values of the 
parameters, and then use the aforementioned map be- 
tween the energy density and the scale factor so that we 
can finally obtain asymptotic expressions for the expan- 
sionary behaviour of the models. Then, we will identify 
the specific late-time behaviour of the models, focusing 
on the existence of future singularities of various types. 
This classification resorts to earlier works by ourselves. 

A separate treatment of the cases = and ^ 
cases is required, so we split the discussion into two 
subsections. 



A. Absent critical density 

In the case of a theory with no critical density, i.e. 
density p* = 0, expressions get considerably simplified: 



-3V^p««/2 
3 



Pit) 



,-3v^(t-to) 



-2/?o 



for Co ^ 0, 



(7) 



for Co = 0. 

Correspondingly, in terms of the expansion factor we get 



a{t) ~ Vk { ^io^o{t-to) 



2/35o 



which provides the following expected results: 

• ^0 < 0: As matter density decreases smoothly, an 
eventual blow up of the corrections to the Fried- 
mann equation is approached. At a finite time 
the scale factor becomes infinite, and the Universe 
experiences a type of singularity which has been 
called "big rip" [l6( (type I in the classification in 

El)- 



4 



Co 


Tipler 


Krolak 


N.O.T. 


(-00,0) 


Strong 


Strong 


I 


[0,oo) 


Non-singular 


Non-singular 


Non-singular 



TABLE II; Singularities in models without critical density 



DGP cosmologies [32|, |33|, |3J, |35|, l36| provide another 
interesting set of examples. If the brane has no tension 
and the bulk is the Minkowski spacetime, one has [13, [131 



(12) 



• $0 > 0: The matter density decreases and the scale 
factor increases smoothly as t grows towards infin- 
ity. This case comprises both quintessence-like be- 
haviors for ^0 G (0,2/3), and non-accelerated evo- 
lutions for ^0 > 2/3. 

• Co = 0: The lowest order term is that of a cos- 
mological constant, and we have to resort to the 
first correction with a positive exponent ^i, which 
leads again to an expression solvable as a Bernouilli 
equation: 



P_ 
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(8) 



In this case 
a{t) ~ ^|e3«i^(*-*«) 



1 hi 
2h^ 
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(9) 



so, this situation represents an exponential expan- 
sion of the Universe, with a corresponding expo- 
nential decrease of matter density, with no future 
singularity at all. 

Therefore, in the case = 0, the modifications consid- 
ered do not produce a qualitative change of behavior to- 
wards the future, except for dramatic modifications pro- 
duced by negative exponents, which lead to a "big rip" 
singularity in the future. 

We close this subsection with several examples which 
fit in case of the general perturbative expression 
of we started from, namely they satisfy 



= hop^° + hip"-^ + 



(10) 



These results are summarized in Table II. 
The first case we consider for illustration is that of the 
power-law Cardassian models [3l| , for which 



= -Gup 1 



P 



(11) 



This expression can be accommodated into (|10l) with 
the following identifications between our parameters and 
those of the original reference: Co = l7Ci="<2/3, 
ho = 87rG/3, hi = (87rG/3)pJj^ > 0. The constant 
Pcard signals the amount of matter energy density p be- 
low which the Cardassian corrections start to dominate 

{Pcard ~ P) 



Here Hq is the value of the Hubble factor today, and 
il^c is the present value of the fractional energy density 
associated with the scale at which the crossover to a cor- 
rections dominated regime occurs. In the perturbative 
formulation required for the discussion we get 



^^(1±1)T ^^'^^''^'^^ 



_ _ (13) 
The self-accelerating branch |33 . |34| arises by taking the 
upper signs, and it is characterized by = and d = 1, 
whereas for the so called normal branch [13, [1H| , which 
arises by taking the lowers signs, one has = 2. 

Finally, we can bring about the reinterpretation of 
Chaplygin-like cosmic evolutions as a modified gravity 
proposal [1^, \^ . In these frameworks one has 



^ + p7(a+l) 



(14) 



with a > and A = (3iJoV(87rG))(i+")(l - , 
and 7 = 1 in the models considered in [39|, whereas 
in the more general framework of [40| 7 is free. Two 
cases are to be distinguished. In the 7 > case, the 
correspondence up to order 7(1 -I- a) in p is given by 
ho = (8G7r/3)Ai/(i+"), = 0, /ii = (87rG/(3(l -I- 
Q!))A~"/(i+"), ^1 = 7(a + 1). But in the 7 < the 
correspondence is rather different, as the identification 
up to order 7 is given by ho = BGtt/S and = 7- 

In all the examples but the last one, according to the 
discussion above, no singular fate of the universe is faced. 
On the contrary, in the last kind of models the singularity 
is of "big rip" type. 



B. Non-trivial critical density 

New features appear for general modifications endowed 
with a non-trivial critical density p*. For this case, we 
assume the matter density has an expansion around the 
critical value = p{to) at a time 

p{t) = + p^(to - ty'^ + P2ito - ty'^ + ■ ■ ■ , (15) 

and from the latter we obtain 

pi_ 

ip* 



a{t) 



1/3 



1 



(16) 



so that the first exponents are the same in both expan- 
sions. 



1, 
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and we may drop the tildes. 
At lowest order we have, 



P ^ 
P 

which upon the requirement of compatibility with Eq. 
^ fixes the first exponent as 



2-6 

The following three cases are to be distinguished: 



(17) 



• 6 < 0: Since < r/i < 1, according to [l5| or Table 
I, these models have a singularity at to with diver- 
gent H (a "big freeze" or singularity type III [11]), 
which is a weak curvature singularity according to 
Tipler [28| , but strong according to Krolak |29|] . 

• 6 € (0,2): In this case 771 > 1, so these models 
could show a weak singularity at tn according to 



15 1 or Table I (sudden singularity [18| or type II in 
13| . or even type IV if ryi > 2, ^ > 1). 



• 6 = 0- The cosmological constant term is domi- 
nant against modifications of the Friedmann equa- 
tion. At first order, we have 



mpi 

p* 



that is, we find a linear behavior for matter density: 

Vi = 1; Pi = ip*\/%)- 

This being so, it turns out we have to expand the 
equation a bit further in order to reveal new quali- 
tative behavior: 
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P 
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ho - ^3\/hop*{to - i)| 



/^-|^(p-p.)«^ + 
2 Vho 



— — [SVhop, 
2?72 ho 



Necessarily, 

?72 = l+6, P2 

and therefore, according to [ll] or Table I 

a{t) = 



(18) 
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(19) 



there is a singularity at to due to the lack of smooth- 
ness of the density and the scale factor. But this 
singularity is weak in both Tipler's [SS] and Krolak' 
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TABLE III: Singularities in models with critical density 



[29| classification, so it does not exert any infinite 
distortion on finite objects going through it and 
cannot, therefore, be considered as a final stage of 
the Universe. It is a sudden singularity or type II 
in [III) for which the scale factor and the density 
remain finite, but H blows up. 

It is worthwhile mentioning that milder singulari- 
ties for which H and also H are finite (type IV in 
[isj ) could, in principle, appear within this frame- 
work, but they would involve choosing 772 > 2, and 
thereby ^ 1, so that the linear term in the den- 
sity in Friedmann equation would be absent. 

Obviously models with analytical expansion, that is, 
natural exponents 6, d,. . . (such as LCDM, for instance) 
do not show future singularities, neither weak nor strong. 

These results are summarized in Table III. 

The normal branch of DGP cosmologies provide a rel- 
evant example for this section. If the bulk on which the 
brane lives is an anti-de Sitter spacetime one has 



H 
i5o 



/SGttp- |A{, 



(20) 



The identification with our perturbed formulation is 
given by = (|Ab| - 3ggr !.J/(8^G), ho = n^^Hl 
^0 = 0, /ii = -4i/oy2^rGa~/(V3) and ^1 = 1/2. This 
sing ularity is a sudden one, also referred to as quiescent 
|38| . or using our terminology, it is a weak extremality 
event. A slight variation leading to a singularity of the 
same sort consists in letting the brane have a negative 
brane tension cr. In this case the bulk can either be the 
Minkowski or the anti-de Sitter spacetime. The above 
expression can be adapted to this variation by simply 
letting p ^ p + u and — > -I- cr. 

Finally, we may consider models arising in loop quan- 
tum cosmology as those in [4111, for which 



SttG 



1-^ 
p* 



but for these models the critical density is relevant for 
the high density regime, imposing a maximum density 
which is reached as the energy density grows. As this 
is the opposite of our working hypothesis (remember we 
demanded p < 0) these models do not quite fit in our 
description here, but could be treated in an analogous 
way, with the corresponding adjustements. 
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III. DISCUSSION 

We here put forward a detailed classification of the 
future behavior of FRW cosmologies in modified grav- 
ity proposals. Departures from the standard description 
of the expansion of the Universe according to Einstein's 
theory have been considered of interest, as they could 
provide an explanation of what is the agent responsible 
for the accelerated expansion of the Universe. 

The main question we pose is what are the characteris- 
tics of the modifications in connection with the presence 
of a singular future behavior of the Universe. As we have 
reflected here, not all the relevant properties of cosmic 
evolution emerge by considering curvature scalars, and 
the deeper insight provided by the consideration of ob- 
servers is needed. 

The spirit of the modified gravity proposals we con- 
sider is to assume the Universe is simply filled with cos- 
mic dust, and no blueshifting component whatsoever is 
considered (unlike when one assumes the current cosmic 
acceleration is due to an exotic fluid or dark energy). Our 
starting point is a perturbative low-energy or infrared ex- 
pansion of the modified Friedmann equation. Two classes 
emerge: those with a critical energy density and those 



without it. We find one has to consider at most the ex- 
ponents of the first two terms of the expansion in order 
to differentiate the possible behaviors, and, more impor- 
tantly, whether the future singularity, if it exists, is weak 
or strong. 

The scheme we propose provides an easy route to con- 
clude the sort of singular behavior present in potential 
new candidates to explain the current acceleration in the 
universe in terms of a modification of gravity. The classi- 
fication we put forward is complementary to others, but 
provide a deeper insight and allow an important further 
degree of refinement. 
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